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Abstract 

It is shown, that the possibihty of a freezing of QCD running cou- 
phng constant at zero in the approach with "forced analyticity" can 
not be in accord with Schwinger-Dyson equation for gluon propagator. 
We propose to add to the analytic expression the well-known infrared 
singular term as well as pole term corresponding to "excited 

gluon" . With this example we formulate the principle of minimality 
of nonperturbative contributions in perturbative (ultraviolet) region, 
which allows us to fix ambiguities in introduction of nonperturbative 
terms and maintain the finiteness of the gluon condensate. As a result 
we obtain estimates of the gluon condensate, which quite agree with 
existing data. The nonzero effective mass of the "excited gluon" leads 
also to some interesting qualitative consequences. 
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The discovery of the asymptotic freedom property |]T| in non-Abehan 
gauge theories turned to be a decisive factor in the formation QCD as the 
strong interaction theory. The negative sign of QCD /3-function = 
Pog^ + ...,Po = — 6o/(167r^), 60 = IIC2/S — 2iVj/3 in the vicinity of zero pro- 
vided the number of active quarks being not too large (for SUc{S) Nf < 16) 
gives couphng constant, which describes quarks and gluons interaction at 
large Euclidean g^, i.e. at small distances, 

^'(^'/^''^)=l-/3o/ln(g7/^^)' 

tending towards zero. Therefore in the deep Euclidean region we are allowed 
to use perturbation theory. In expression (0), which takes into account the 
main logarithms, /i is a normalization point. An account of the next correc- 
tions does not change asymptotic behaviour (|I]) for ^ 00. By introducing 
dimensional constant A^ = /i^ exp(— 47r/(6oas)), c^s = 5'^/47r, we turn from 
explicitly renormalization invariant expression (|l]) to the following formula 

It is reasonable to estimate parameter A in approximate expression (|^) to 
be around of few hundreds MeV. With decreasing effective constant (|^) 
increases, that may indicate a tendency of unlimited growth of the interaction 
at large distances, leading to a confinement of coloured objects. However, at 
g2 — ^2 expression (|^) the pole is present, which is nonphysical at least 
due to failing of the perturbation theory, starting from which formula @) has 
been obtained. 

In recent work a solution of the problem of ghost pole was proposed 
with imposing of a condition of analyticity in q^. The idea of "forced ana- 
lyticity" goes back to works [|, ^ of the late fifties, which were dedicated 
to the problem of Landau-Pomeranchuk pole in QED. Using for oisiq^) 
a spectral representation without subtractions, the following expression for 
the running coupling constant was obtained in paper [0] 



A2 



(3) 



In(g2/A2) A2 - q^ 

This expression has the asymptotic freedom property and its analyticity in 
the infrared region is due to nonperturbative contributions. It does not con- 
tain any additional parameter and has finite limit at zero, q;^^-*(0) = 47r/6o — 
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1.40 (freezing of the coupling constant), which depends only on symmetry 
factors. This limit turns out to be stable with respect to higher orders cor- 
rections. 

As it is noted in work a procedure of summation of leading loga- 
rithmic terms is not defined uniquely. A partial fixation of this ambiguity 
in QED is realized by using of a method of summation of the perturbation 
theory series under the sign of the spectral integral of the Kallen-Lehmann 
representation. Nevertheless after such summation a functional ambiguity 
remains, which from one side does not violate correct analytic properties 
of Green functions in a complex plane of a corresponding invariant variable 
and from the other side contains nonanalytic dependence on constant g^. In 
work 1^ while investigating the photon propagator in QED it was shown, that 
ambiguities in summation procedure of the diagram series can be removed 
provided one demands not only the validity of spectral representation, but 
also the fulfillment of equations of motion. 

In the present note we consider a problem of consistency of the constant 
behaviour of the effective charge in the infrared region with Schwinger-Dyson 
(SD) equation for a gluon propagator. Further we include into consideration 
nonperturbative terms, the singular in the infrared region term ~ in 
particular, the necessity of the renormalization invariance being taken into 
account. Then we discuss possibilities of an adjustment of demands of con- 
finement, asymptotic freedom, analyticity, an accordance with the pertur- 
bation theory and correspondence with estimates of the gluon condensate 
value. 

To study the problem of a possibility of a constant behaviour of the run- 
ning constant in the infrared region let us consider the integral SD equation 
for the gluon propagator in ghost-free axial gauge |0] A'p]^^ = 0, — gauge 
vector, rf ^ In this gauge the effective charge is directly connected with 
the gluon propagator and Slavnov- Taylor identities have the most simple 
form. The important preference of the axial gauge consist in the possibility 
to exclude the term from SD equation, which contains the full four-gluon 
vertex by means of contraction of the equation with tensor rjf^rjy/rf. 

In what follows we shall work in the Euclidean momentum space, where 
smallness of the momentum squared is immediately connected with smallness 
of its components. The equation to be considered has the form: 



[^,j(p)-^roU(p)] 
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„2. 4-n „ 

= —§£yr j d'-krf^xM ^ - P)DxAk)DAp -k)x 

xT3^Su{k,p- k,-p), (4) 

where Il^^{p) is the one-loop part of the polarization operator, D^y{p) is 
the propagator, V„s^{k^p — k, —p) is the one-particle irreducible three-gluon 
vertex function, T^^lp{p, —k, k — p) is the free three-gluon vertex. 

We suppose the approximation D^^^p) = Z{p^)D((j)^y{p) to be appropri- 
ate to study the infrared region. Let us divide the momentum integration 
domain in expression (Q) in two parts: k"^ < and fc^ > A^, where A is 
sufficiently small, but finite. Then domain k"^ > A^ in the case of absence 
of kinematic singularities in three-gluon vertex gives a contribution, which is 
regular in p'^ for ^ 0, and in domain k'^ < A^ full Green functions can be 
approximated by free ones up to constant factors according to an assumption 
of running constant be frozen at zero. Then one can write 

nA«-(P)^ = y ^ kT',>,pip,-k,k-p)x 

X - ^)r£L(fc,P - k, -p)v,vu/v' + Q{p'; y, K n). (5) 

Here y = {pviY /p'^rf is the gauge parameter. The integration in formula (|^) 
can be extended up to all the domain of momentum, that results in a change 
of the regular in p^ contribution Q. Thus one has 

n,.(p)^ = ^(0)n2(p)^ + Q{p'; y. n), (6) 

where n|/J(p) is the one-loop perturbation theory contribution to the polar- 
ization operator. This contribution is calculated in paper and has rather 
complicated structure. Let us present the expression for the leading terms of 
convolution (BI) at ?/ — 0. We have 



mp)'-f = cp' 




70 



40 , ^, 2^ ^ 



(7) 
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Here C = (7^C2/327r^, 7 is the Euler constant. From expression (|^ we see, 
that singularity at y = is smooth and the hmit at y = does exist. Term 
~ 1/e (n = 4 + 2e) as well as constant ones could be absorbed into function 
Q while the logarithm of the momentum squared necessarily persists. The 
equation for function Z{p'^) takes the form 
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n] 



(8) 



does not agree with the 
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We see, that behaviour ZijP') ~ Z{Q) 7^ for 
SD equation. 

This conclusion stimulate us to look for the possibilities different from 
the assumption on the finiteness of the coupling constant at zero. Recently 
the possibility of the soft singular power infrared behaviour of the gluon 
propagator was discussed |T^, D{q) ~ {<f')~^^ — ^ 0, where /? is a small 
positive non-integer number. In Ref. [|ll|] the consistency of such behaviour 
with Eq. (Q) was studied. A characteristic equation for the exponent (5 was 
obtained and this equation was shown not to have solutions in the region 
< /? < 1. The authors of Ref. |T2| also come to the conclusion on the 
inconsistency of the soft singular infrared behaviour of the gluon propagator. 
The case of possible interference of power terms was studied in Ref. [0] and it 
was shown that in a rather wide interval — 1 < /3 < 3 of the non- integer values 
of the exponent the characteristic equation has no solutions. At present the 
more singular, in comparison with free case, infrared behaviour of the form 
D{q) ^ M^/(g^)^, — > seems to be the most justified 0, |15|, |T6|. The 
physical consequences of such enhancement of zero modes are discussed in 
the reviews [17, 18 1. Bearing in mind the remarks stated above let us consider 



the following expression for the running coupling: 



Air 



+ 



A' 



+ c- 



(9) 



Ing2/A2 ' A^ — 

Let us represent this expression in explicitly renormalization invariant form. 
It can be done without solving the differential renormalization group equa- 
tions. In this order we write as(g^) = ^^(g^//i^, fi'^)/47r and use the normal- 
ization condition g'^{l,g'^) = g^. Then we obtain the equation for wanted 
dependence of the parameter A^ on g"^ and /i^: 



g^An 



+ 



A' 



In /i2/A2 A"^ - jj 



+ c- 



A^ 
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From dimensional reasons = fi'^exp{—(f{x)}, where x = bog'^/lGir'^ = 
boas/ 4:71, and for function (p{x) we obtain the equation: 

^ ~ Lp{x) 1 - e'^(^) '^^ 

The solution of this equation at c > is monotone decreasing function f{x), 
which has the behaviour (p{x) ~ 1/x at a: — >■ and if{x) ~ — ln(a;/c) at 
X — s> +00. The relation obtained ensures the renormalization invariance of 
as(g^). At low we obtain A^ = /x^ exp{— 47r/(6otts)}, which indicates the 
essentially nonperturbative character of both last terms of the Eq. @ and 
this terms are absent in the perturbation theory. With given value of the 
QCD scale parameter A the parameter c can be fixed by the string tension 
K or the Regge slope a' = l/{2TrK) assuming the linear confinement V{r) ~ 
Kr = o?r at r — > oo. We define the potential V{r) of static qq interaction 
1^, ^ by means of three-dimensional Fourier transform of as{(f)/(f with 



the contributions of only one dressed gluon exchange taken into account. 
This gives the following relation 

cA^ = (36o/87r)a^ = (6o/167r^)/M^ (10) 

At large q^ from Eq. (Rl) one obtains 



«sW ) = -7- 



+ (c-i)^-74^ + 0((gV) 



Ing^/A^ q^ [q^ 



From Eq. (0) it is seen that in the ultraviolet region the nonperturbative 
contributions decrease more rapid then all renormalization group improved 
perturbation theory corrections. The value c = 1 corresponds to maximal 
suppression of nonperturbative contributions in the ultraviolet region. Ac- 
cepting this condition one obtains the connection of the QCD scale parameter 
A and string tension /t = of the form A^ = 36o/^/87r . Taking a ~ 0.42 GeV 
one obtains for A reasonable estimation, A ~ 0.434 GeV {Bq = 9 in the case 
of 3 light flavours). 

Considering the nonperturbative contributions the following arguments 
can be expressed. One knows QCD to be renormalizable in the perturbation 
theory and, as usually, the renormalization procedure can be developed to 
remove the divergences in all orders. However, what about the nonpertur- 
bative contributions? If they bring in the additional divergences then the 
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problem of renormalization turns out to be unsolved. The situation when 
nonperturbative contributions do not violate the perturbative renormaliza- 
tion properties seems to be more attractive. It take place if the nonperturba- 
tive contributions decrease at momentum infinity sufficiently fast and do not 
introduce the divergences in observables. So, it is natural to demand their 
fastest of possible decrease at large momenta. An application of the principle 
of minimality of nonperturbative contributions in the ultraviolet region will 
be shown further with taking as an example the important physical quantity, 
namely, the gluon condensate, K = < a^/vr : G^^, :> . According to the 
definition (see e.g., Jl^) up to the quadratic approximation in the gluon fields 
one has after Wick rotation 

^ = V / (2^ (^^('') - ^'"('')) = vJo ^^"^y^y^y^ ^'^^ 

where a^'^ is nonperturbative part of the running coupling constant. In our 
case the two last terms of Eq. (||) should be taken. By substituting this terms 
in Eq. ([T^) one can see the logarithmic divergences of the integral at infinity 
and at finite point k"^ = h?. 

The acceptance of the cancellation mechanism for the nonphysical per- 
turbation theory singularities (|) by the nonperturbative contributions leads 
to the necessity of supplementary definition of the integral (0) near point 
k'^ = h?. This problem can be reformulated as a problem of dividing of per- 
turbative and nonperturbative contributions in resulting in introduction 
of some parameter /cq = 1 -i- 2 GeV. This provides absence of the pole at 
/c^ = A^ in both perturbative and nonperturbative parts. The divergence of 
the integral (|I2]) at infinity stimulate the further modification of the running 
coupling constant. Going over from Eq. (^ to Eq. (^ the isolated singularity 
was introduced. In this case the singularity corresponding to the unitary cut 
was not changed and in accordance with the approach of Refs. ^ |, || is de- 
termined by perturbation theory. Following to this logic let us consider the 
expression for cxs with one more isolated singularity in the time-like region. 
The tachion singularity in the space-like region, of cause, is prohibited. 

The principle of minimality of nonperturbative contributions in ultra- 
violet region then leads to the following unique expression for the running 
coupling constant 

_ (^2) _ 47r/ 1 ^ A^ + — + ~ ^)^' 
6o \ln{q^ / K^) Is? — q^ q"^ + rrig 
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with fixed residue and mass parameter m^, 

mJ = AV(c-l), (14) 

for tlie newly introduced term. Tlie expression (0) can be represented in ex- 
plicitly renormalization invariant form in a similar way to the expression (|^). 
Nonperturbative contributions in Eq. (p!3D decrease at infinity as l/g^, the 
integral in Eq. ([l^ ) converges and we can obtain 

K = -1- {ln{c - 1) + kl/K^ + ln{kl/K^ - 1)} . (15) 

Phenomenology gives the positive value of the gluon condensate K in the 
interval (0.32 Gel/)^ - (0.38 Gel/)^ [0, 0. As an example we take values 



= 1.2 ^ 1.3 GeV. If one regards the string tension parameter to be given, 
then from Eqs. (|T^, ( [T5| ) and ( PD] ) one has the dependencies of all the values 



under consideration on the parameter c, which are presented in the Table I. 
Note that values c = 1.063, A = 422 Me V, rUg = 1.682 Gel/, ko = 1.265 



GeV corresponds to the conventional value of the gluon condensate |21| K = 
{0.33 GeV)'^. Certainly these results should be considered as tentative, but 
nevertheless they seems encouraging. 

It is seen from Eq. (|13D that the pole singularities are situated at two 
points = and = — "n^g- It corresponds to two effective gluon masses, 
and TUg. Therefore the physical meaning of the parameter is not the 
constituent gluon mass but rather the mass of the exited state of the gluon. It 
is essential that the residue at is very small, so the states with the exited 
gluons should be quite narrow in contrast to the spectrum of the coupled 
massless gluons. 

The qualitative picture of the glueball states corresponding to the running 
coupling constant (|T3|) with nig ~ 1.7 GeV could be the following: 

1) The states gg — continuous spectrum and very wide resonances are 
probable; 

2) The states g g' — resonances with probable mass interval 1500 - 1800 
MeV and with width suppression factor (1 — c); 

3) The narrow states g' g' — resonances with possible masses 3000 - 3600 
MeV and with width suppression factor (1 — c)^. 

Note that in the region 2) there are the glueball candidates. The region 3) 
is insufficiently investigated, some indications in favour of the narrow states 
are showing up (see e.g., [p3l). 
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Table 1: Parameters of the running coupling constant (13) and gluon con- 
densate as functions of parameter c. 



c 


A, GeV 


nig, GeV 


GeV 


GeV 


K^^ GeV 






ko = 1.2 GeV 


ko = 1.25 GeV 


ko = l.3 GeV 


1.01 


0.433 


4.332 


0.298 


0.309 


0.318 


1.02 


0.431 


3.048 


0.307 


0.317 


0.326 


1.03 


0.429 


2.476 


0.312 


0.321 


0.330 


1.04 


0.427 


2.134 


0.315 


0.324 


0.332 


1.05 


0.425 


1.900 


0.317 


0.326 


0.334 


1.06 


0.423 


1.726 


0.319 


0.327 


0.335 


1.07 


0.421 


1.591 


0.320 


0.328 


0.336 


1.08 


0.419 


1.481 


0.321 


0.329 


0.337 


1.10 


0.415 


1.313 


0.322 


0.330 


0.337 


1.12 


0.411 


1.187 


0.323 


0.330 


0.337 


1.16 


0.404 


1.010 


0.323 


0.330 


0.337 


1.20 


0.397 


0.889 


0.322 


0.329 


0.336 


1.24 


0.391 


0.798 


0.321 


0.328 


0.335 


1.30 


0.382 


0.697 


0.319 


0.326 


0.332 
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